lOSR Journal of Engineering (lOSRJEN) 

ISSN (e): 2250-3021, ISSN (p): 2278-8719 

Vol. 05, Issue 02 (February. 2015), IIV4II PP 01-12 



www.iosrjen.org 



PH/PH/1 Bulk Arrival and Bulk Service Queue with Randomly 

Varying Environment 

Ramshankar.R\ Rama.G^, Sandhya.R^, Sundar.V^, Ramanarayanan.R^ 

'independent Researcher MS9EC0, University of Massachusetts, Amherst, MA, USA 
^Independent Researcher B. Tech, Vellore Institute of Technology, Vellore, India 
^Independent Researcher MSPM, School of Business, George Washington University, Washington .D.C, USA 
"^Senior Testing Engineer, ANSYS Inc., 2600, Drive, Canonsburg, PA 153 17, USA 
^Professor of Mathematics, (Retired), Vel Tech University, Chennai, INDIA. 



Abstract: This paper studies two stochastic bulk arrival and bulk service PH/PH/1 queue Models (A) and (B) 
with randomly varying k* distinct environments. The arrival and service distributions are {ai,Ti) and (Pi,Si) in 
the environment i for 1 <i<k* respectively. Whenever the environment changes from i to j the arrival PH and 
service PH distributions change from the i version to the j version with the exception of the first remaining 
arrival time and first remaining service time have stationary PH distributions of the j version which is known as 
equilibrium PH distribution for 1 < i, j < k* and on completion of the same the arrival and service distributions 
become initial versions {aj,Tj) and (J3j,Sj^. The queue system has infinite storing capacity and the state space 
is identified as five dimensional one to apply Neuts ' matrix methods. The arrivals and the services occur 
whenever absorptions occur in the corresponding PH distributions. The sizes of the arrivals and the services are 
finite valued discrete random variables with distinct distributions with respect to environments and with respect 
to PH phases from which the absorptions occur. Matrix partitioning method is used to study the models. In 
Model (A) the maximum of the arrival sizes is greater than the maximum of the service sizes and the 
infinitesimal generator is partitioned mostly as blocks of the sum of the products of PH arrival and PH service 
phases in the various environments times the maximum of the arrival sizes for analysis. In Model (B) the 
maximum of the arrival sizes is less than the maximum of the service sizes. The generator is partitioned mostly 
using blocks of the same sum-product of phases times the maximum of the service sizes. Block circulant matrix 
structure is noticed in the basic system generators. The stationary queue length probabilities, its expected 
values, its variances and probabilities of empty levels are derived for the two models using matrix methods. 
Numerical examples are presented for illustration. 

Keywords: Bulk Arrivals, Bulk Service, Block Circulant Matrix, Neuts Matrix Methods, Phase Type 
Distribution, Stationary PH Distribution. 



1. INTRODUCTION 

In this paper two bulk arrival and bulk service PH/PH/1 queues with random environment have been 
studied using matrix geometric methods. Numerical studies on matrix methods are presented by Bini, Latouche 
and Meini [1]. Multi server model has been of interest in Chakravarthy and Neuts [2]. Birth and death model 
has been analyzed by Gaver, Jacobs and Latouche [3]. Analytic methods are focused in Latouche and 
Ramaswami [4] and for matrix geometric methods one may refer Neuts [5]. For M/M/1 bulk queues with 
random environment models one may refer Rama Ganesan, Ramshankar and Ramanarayanan [6] and M/M/C 
bulk queues with random environment models are of interest in Sandhya, Sundar, Rama, Ramshankar and 
Ramanarayanan [7]. PH/PH/1 bulk queues without variation of environments have been treated by Ramshankar, 
Rama Ganesan and Ramanarayanan [8]. The models considered here are general compared to existing models. 
Random number of arrivals and random number of services are considered at a time. Fixed numbers of 
customers are cleared by a service in the models of Neuts and Nadarajan [9]. In real life situations when a 
machine manufactures a fixed number of products in every production schedule, the defective items are rejected 
in all production lots', making the production lot is only of random size and not a fixed one always. Situations 
of random bulk services are seen often in software based industries where finished software projects waiting for 
marketing are sold in bulk sizes when there is economic boom and the business may be insignificant when there 
is economic recession. In industrial productions, bulk types are very common. Manufactured products arrive in 
various bulk sizes for sale in markets and the products are sold in various bulk sizes depending on market 
requirements. Noam Paz and Uri Yechali [10] have studied M/M/1 queue with disaster. Usually bulk arrival 
models have M/G/1 upper-Heisenberg block matrix structure. The decomposition of a Toeplitz sub matrix of the 
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infinitesimal generator is required to find the stationary probability vector as done in William J. Stewart [11] 
and matrix geometric structures have not been noted. In such models the recurrence relation method to find the 
stationary probabilities is stopped at a certain level in most general cases using a terminating analysis very well 
explained by Qi-Ming He [12] and this stopping limitation of terminating method converts an infinite arrival 
system to a finite arrival one. In special cases generating function has been identified by Rama and 
Ramanarayanan [13]. However the division modulo partitioning of the infinitesimal generator along with 
environment and PH phases used in the paper is presenting matrix geometric solution for finite sized arrivals 
and services models. The M/PH/1 and PH/M/C queues with random environments have been studied by Usha 
[14] and [15] without bulk arrivals and bulk services. It has been noticed by Usha [14, 15] that when the 
environment changes the remaining arrival and service times are to be completed in the new environment. The 
residual arrival time and the residual service time distributions in the new environment are to be considered in 
the new environment at an arbitrary epoch since the spent arrival time and the spent service time have been in 
the previous environment with distinct sizes of PH phase. Further new arrival time and new service time from 
the start using initial PH distributions of the new environment cannot be considered since the arrival and the 
service have been partly completed in the previous environment indicating the stationary versions of the arrival 
and service distributions in the new environments are to be used for the completions of the residual arrival and 
service times in the new environment and on completion of the same the next arrival and service onwards they 
have initial versions of the PH distributions of the new environment. The stationary version of the distribution 
for residual time has been well explained in Qi-Ming He [12] where it is named as equilibrium PH distribution. 
Randomly varying environment PH/PH/1 queue models with bulk arrival and bulk service have not been treated 
so far at any depth. In this paper the partitioning of the matrix is carried out in a way that the stationary 
probability vector exhibits a matrix geometric structure for PH/PH/1 bulk queues with random environment 
where the arrivals and service sizes are finite. Two models (A) and (B) on PH/PH/1 bulk queue systems with 
infinite storage space for customers are studied using the block partitioning method. Model (A) presents the case 
when M, the maximum of the arrival sizes is bigger than N, the maximum of the service sizes. In Model (B), its 
dual case N is bigger than M, is treated. In general in Queue models, the state space of the system has the first 
co-ordinate indicating the number of customers in the system but here the customers in the system are grouped 
and considered as members of blocks of sizes of the maximum for finding the rate matrix. Using the maximum 
of the bulk arrival size or the maximum of the bulk service size and grouping the customers as members of 
blocks in addition to coordinates of the arrival and service phases for the partitioning the infinitesimal generator 
is a new approach in this area. The matrices appearing as the basic system generators in these two models due to 
block partitioned structure are seen as block circulant matrices. The paper is organized in the following manner. 
In sections II and III the stationary probability of the number of customers waiting for service, the expectation 
and the variance and the probability of empty queue are derived for these Models (A) and (B). In section IV 
numerical cases are presented to illustrate them. 

II.MODEL (A): MAXIMUM ARRIVAL SIZE M > MAXIMUM SERVICE SIZE N 
2.1 Assumptions 

(i) There are k* environments. The environment changes as per changes in a continuous time Markov chain with 
infinitesimal generator of order k* with stationary probability vector n'. 

(ii) In the environment i for 1 < i < k*, the time between consecutive epochs of bulk arrivals of customers has 
phase type distribution ( cfj, 7;) where Tjis a matrix of order fe; with absorbing rate 7'; = —Tie to the absorbing 
state fej-Hl from where the arrival process moves instantaneously to a starting state as per the starting vector cf; 

= {ai l, '^1,2-' '^i.ki) ^yli '^i.j ~ 1- L^*- ^ invariant probability vector of the generator matrix 
(Ti + T'iai). 

(iii) In the environment i for 1 < i < k*,when the absorption occurs in the PH arrival process due to transition 
from a state j to state fej +1, jX number of customers arrive with probabilities P (jj= n) = jPn for 1 < n < jM 

and Zin=i jp„=l where jM is the maximum arrival size for PH phase j where 1< j < fej . 

(iv) In the environment i for 1 < i < k*, the time between consecutive epochs of bulk services of customers has 
phase type distribution (Pi, Si) where 5; is a matrix of order k'l with absorbing rate S'l = —SiC to the 
absorbing state k' i+\ from where the service process moves instantaneously to a starting state as per the staring 
vector Pi = {Pix>Pi,2- > ■■■> Pi.k'i) ^y=i ~ ^- ^ invariant probability vector of the generator 
matrix (5; -I- 5'^ 

(v) In the environment i fori < i < k*, customers of bulk size jtp are served at epoch when the absorption occurs 

due to a transition from state j to state k'i+1, with probabilities P(jT/; = n) = jq„ for 1< n < jN and Y/j^^ jQn = 
1, when more than jN customers are waiting for service where jN is the maximum service size for PH phase j 
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where 1 < j < fc'j . When n customers n < jiV are waiting for service, then n' customers are served with 
probability jq^- for 1< n' < n-1 and n customers are served with probability z/.^^jqn for PH phase j where 

i<j<fcV 

(vi) When the environment changes from i to j for 1 < i, j < k*, the arrival and service distributions in the new 
environment j are the stationary (equilibrium )versions of arrival time and service time distributions in the new 
environment, namely, {q)j , Tj ) and {(pj , Sj ) respectively for the completions of the residual arrival and service 
times and on completion of the same the next arrival and service onwards they have initial versions of the PH 
distributions of the new environment namely (a, , Tj ) and (Pj , Sj ) respectively. 

(vii) The maximum arrival size M= maXi^^^ maxj^ <, jM is greater than the maximum service size 
N= maxi<i<fc, maxi <,■ jiV . 

2.2.Analysis 

The state of the system of the continuous time Markov chain X (t) under consideration is presented as follows. 
X(t)= {(0, i,j) : for 1 <i<k* 1 <j </£;)} U {(0, k, i,j,j') ; for 1 <k<M-l; 1 < i < k*; 1 < j < fe;; 1 <}<k'i} 
U {(n,k,i,j,j'): forO<k<M-l; 1 < i < k*; 1 < j < /c^; 1 <j<fe'i andn> 1}. (1) 
The chain is in the state (0, i, j) when the number of customers in the queue is 0, the environment state is i for 1 
< i < k*and the arrival phase is j for 1 < j < fcj. The chain is in the state (0, k, i, j, j') when the number of 
customers is k for 1 < k < M-1, the environment state is i for 1 < i < k*, the arrival phase is j for 1 < j < fe; and 
the service phase is j' for 1 < j' < fe';. The chain is in the state (n, k, i, j, j') when the number of customers in the 
queue is n M + k, for 0 < k < M-1 and 1 < n < co, the environment state is i for 1 < i < k*, the arrival phase is j 
for 1 < j < fej and the service phase is j' for 1 < j' < fe';. When the number of customers waiting in the system is 
r, then r is identified with (n, k) where r on division by M gives n as the quotient and k as the remainder. Let the 
survivor probabilities of arrivals]/ and of services jXp be respectively P(jx>m)= jP^ =1-I]^i jPn^for 
1 <m<]M-l and l<j<fei (2) 
P(lxp>m)=lQ^=l-'Z.jr=ijqn, for 1 <m<;iV-land l<j <fe'i (3) 
with jPq = 1, for all j, 1< j < fe; and jQq = 1 for all j , 1< j < k'l for the environment state i for 1 < i < k* . 
The chain X (t) describing model has the infinitesimal generator of infinite order which can be presented in 
block partitioned form given below. 
rBi Bo 0 0 . . . ■ 



Qa-- 



B2 Aq 0 

0 i4i i4o 0 

0 0 A2 Ai Ao 0 

0 0 0 A2 Ai Ao 



(4) 



In (4) the states of the matrices are listed lexicographically as 0, 1, 2, 3, .... For partition purpose the zero states 
in the first two sets given in (1) are combined. The vector 0 is of type 1 x [I]f=i fe; + (M — 1) Xf=i fe; fe';] and is 

0=((0,1,1),(0,1,2),(0,1,3)...(0,1,/Ci),(0,2,l),(0,2,2),(0,2,3)...(0,2,fe2), (0,k*,l),(0,k*,2),(0,k*,3). . .(0,k*,fefc J, 

(0,1, 1,1,1),(0,1, 1,1,2).... (0,1, l,l,fe'i),(0,l,l,2,l),(0,l, 1,2,2).... (0,1, l,2,fc'i),(0,l, 1,3,1).... (0,1, l,3,fe'i) (0,1,1, 

fci, l)...(0,l,l,fei,fe'i),(0,l,2,l,l),(0,l,2,l,2)... .(0,1,2,1, fc'2),(0,l,2,2,l),(0,l,2,2,2)....(0,l,2,2,fe'2),(0,l,2,3,l)....(0 
,l,2,3,/c'2)....(0,l,2,fe2,l)...(0,l,2,fe2,fc'2),(0,l,3,l,l)...(0,l,3,fe3,fe'3)...(0,l,kM,l),...,(0,l,k*,fcfc.,fe'fc.),(0,2,l,l 
,l),(0,2,l,l,2),...,(0,2,k*,fcfe.,fc'fcJ,(0,3,l,l,l)...(0,3,k*,fefc.fcfcJ,(0,4,l,l,l)...(0,4,k*,fcfc./cV.)...(0,M-l,l,l,l)... 
(0,M-l,k*,fej.j, /c'j.j)) and the vector n is of type lx[Ml]f=ifeife'i] and is given in a similar manner as follows 
nKn,0,l,l,l),(n,0,l,l,2)....(n,0,l,l,fe'l),(n,0,l,2,l),...(n,0,l,2,fe'l),...(n,0,k^l,l),...(n,0,k^fefc.,feVJ,(n,l,l,l,l) 

....(n,l,k*,fcfc.,fc'fcj,(n,2,l,l,l)....(n,2,k*,fcfc.,fe'fc.) (n,M-l,l,l,l),(n,M-l, 1,1,2) (n,M-l,k*,fcfc., fe'^.)). 

The matricesB]^ and have negative diagonal elements, they are of orders [I]f=i fe; + (M — 1) Xf=i fe; ^'J and 
[M Yj'i=i k-i k' i] respectively and their off diagonal are non -negative. 

The matrices Aq andA2 have nonnegative elements and are of order [ M Xf=i fe; k'l ] . The matrices Bq and B2 
have non-negative elements and are of types [Xf=i fe; + (M — 1) Xf=i fe; fe'J x [ MXf=ifeife'i ] and 
[ M Xf=i k^ k'l ] X Ef=i fe; + (M - 1) Yn=i ki fe';]. Component mati'ices of A^ and Sjfor i=0,l,2 ai-e defined 
below. Let © and ® denote the Kronecker sum and Kronecker product. 

Let 2;=r, 05, + diag {{Q^\{) = (7^(8)4-.) + {h,®S{) + diag {{Q^\{) for 1 < i < k* (5) 
where I indicates the identity matrices of orders given in the suffixes, Qi is of order fejfe'j and the last term is a 
diagonal matrix of order fejfe'j . Considering the change of environment switches on stationary distributions in 
PH arrival time and PH service time in the new environment, the following matrix Q of order YJi=iki k'iis 
defined which is concerned with change of environment during arrival time and service time. 
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rQ'i 


^1,2 


^1,3 ■ 


■ ^l,k* 




^2.1 


Q'2 


^2.3 • 


■ ^2,k* 




^3.1 


^3.2 


Q'3 ■ 


■ ^-i.k* 




-^k*,l 


^k*,2 


■^fc*,3 ■ 


■ Q'k.- 



(6) 



where Hij is a rectangular matrix of type fejfe'j x kj fey whose all rows are equal to {Qi)ij {(Pj ® (f>j) for i i^j , 
1 < i, j < k*. The arrival rate of n customers for 1< n < jM corresponding to absorption to state fej+1 from the 
arrival PH phase j for 1< j < fej, in the environment i for 1< i < k* is given by the component j of the column 
vector of type k, xl where r;,„= ( (T\\(ipJ , {T\UiVn) , (T'ih(iPn) , ■■■■(T'dk^kJ^n) )'; (7) 
the service rate of n customers for 1< n < jN corresponding to absorption to state fe'j+1 from the service PH 
phase j for 1< j < fe'j, in the environment i for 1< i < k* is given by the component j of the column vector 5; „ 
of type k\ xlwhere 5;,„= ( (S'Miq,) , {S'Miqn) , (S'tUkn) , ■■■■iS't)k'Sk'icin)y ■ (8) 
[Tin 0^1® Ik' ^ 0 0 



Let A„ 



T2,n «2 ® Ik'2 

0 



0 

Txn «3 ® Ik 



for 1 < n < M (9) 



0 0 0 - T,,,nak.®Ik',. 

In (9) An is a square matrix of order J]f=i fc; k\; Tj „ aj ® /fc' Js a square matrix of order kj k'j for 1 < j < k* and 
0 appearing as (i, j) component of (9) is a block zero rectangular matrix of type fejfe'; x fey fey . 
[I,,® sin Pi 0 0 ■•■ 0 

0 Ik2®S'2,nP2 0 - 0 

Let[/„= 0 0 43®4n/?3 - 0 



for 1 <n<N (10) 



0 0 0 ••■ 4,. ® /?, J 

In (10) Un is a square matrix of order Xi=i fe; fe'ii Ikj ® n Pj^^ ^ square matrix of order fey fey for 1 < j < k* 
and 0 appearing as (i, j) component of (10) is a block zero rectangular matrix of type fej fe'j x fey fey .The matrix Ai 
for i = 0,1,2 are as follows. 



Am 0 

^M-l 

^M-2 ^M-1 
^M-2 



A3 



A4 
A3 
A2 



0 
0 
0 
0 

Am 0 
Am-1 Am 
Am-2 Am-I 



0 
0 
0 
0 

0 
0 

Am 



(11) 



A, = 



-0 ■ 


■ 0 




Un-1 ■ 


■ ^2 


Ui ] 


0 ■ 


■ 0 


0 


Un ■ 




U2 


0 ■ 


■ 0 


0 


0 


■ y« 


Un-1 


0 ■ 


• 0 


0 


0 


• 0 


Un 


0 ■ 


■ 0 


0 


0 


■ 0 


0 


-0 ■ 


■ 0 


0 


0 


■ 0 


0 - 



(12) 



-n 


Ai 


A2 ■ 


■ Am-n- 


-2 


Am-n-1 


Am-n 


Am-2 


Am-1 


Ui 


n 


Ai ■ 


■ Am-n- 


-3 


Am-N-2 


Am-n-1 " 


Am-3 


Am-2 


f/2 


Ui 


n 


■ Am-n- 


-4 


Am-N-3 


Am-n-2 " 


Am-4 


Am-3 


Un 


Un-1 


Un-2 ■ 


n 




Ai 


A2 


■ Am-N-2 


Am-n-1 


0 


Un 


Un-1 ■ 


Ui 




n 


Ai 


■ Am-N-3 


Am-n-2 


0 


0 


Un ■ 


U2 




Ui 


n 


■ Am-n-^ 


Am-N-3 


0 


0 


0 


■ u^ 




Un-1 


Un-2 ■■ 


n 


Ai 


. 0 


0 


0 


0 




Un 


Un-1 ■■ 


Ui 


n 



For defining the matrices Bj for i = 0,1,2 the following component mafrices are required 



A' = 



Ti,„ O ft 0 

0 T2,n «2 ® P2 

0 0 



0 
0 

T3.n «3 ® /?3 



0 



0 



0 
0 
0 

"^fc+.n ^k* 



for 1 < n < M 



(13) 



(14) 



A'n is a rectangular matrix of type ( Xi=i ki)x Hj^i (fej fe'j) for 1 < n < M ; Tj „ ctj (S) ft is a rectangular matrix 
of order fcjXfejfe'j and 0 appearing as (i, j) component of (14) is a block zero rectangular matrix of type fej x fey fe'y 
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for 1 < i, j <k*. Let 7% = h,®{ {S'^^Qn) • (5' JzCiQn) , (5' JsGQn) (fc'!'?n))' 

for 1 < n < N -1 is a matrix of type fejfe'j x fe; for 1 < i < k*and let 
0 0 ■•■ 0 



K, = 



0 



0 



for 1 < n < N. 



(15) 



0 0 0 - n..J 

This is a rectangular matrix of type ( 1];=!* kik'i) x (l]f=i fe;) and 0 appearing in the (i, j) component is a 
rectangular 0 matrix of type fcjfc'j x kj for 1 < i, j < k*. 



LetU = 



fci ® 5i 
0 



0 



k2 



(16) 



0 0 0 - 4^0 5fc. 

In (16), U is a rectangular matrix of type ( I]i=i* fci^'i) ^ (lif=i ''^i) ™d 0 appearing in the (i, j) component is a 
rectangular 0 matrix of type fejfe'j x kj for 1 < i, j < k*. 4. (8)5; is a rectangular matrix of type fejfe'j x fe; for 1 < 
i < k*. The matrix Bq is same as that of Aq when /l^ in the first row of v4Qis replaced byyl'^. The matrix B^ is 
given below. The matrix B2 is same as that of v42when the first block column with 0 is considered as Xf=i fe; 
columns block instead of Xf=i feife'jColumns block of ^42- To write Bj^the block for 0 is to be considered which 
has queue length, L= 0, 1,2.. .M-1. When L = 0 there is only arrival process and no service process. The change 
in environment from i to j switches on stationary PH (equilibrium PH) distribution in the new environment j 
whenever it occurs for 1 < i 7^: j, < k*. When an arrival occurs and queue length becomes L in the environment i 
both the arrival time and the service time start with starting probability vector ct; and respectively for 1 < i < 
k*. In the 0 when L =1,2, . . .M-1 all the processes arrival, service and environment are active as in other blocks 
n for n > 0. Considering the change of environment switches on the stationary (equilibrium) distribution in PH 
arrival time in the new environment when the queue is empty, the following matrix O' of order Xf=i fejis defined 
which is concerned with change of environment during arrival time. 
T\ n'i2 ^'i.-i ■'■ ^'l,k* 



(17) 



^'2.1 
^'3.1 



r2 
^'3.2 



^'2.3 



n' 



k*,l 



n' 



k*,2 



n' 



kt,3 



^' 2,k* 
^\,k* 



Here T'i= 7; + diag{Qi)i i and Q' q is a rectangular matrix of type fc; x kj whose all rows are equal to 
(Qi)i,j % presenting the rates of changing to phases in the new environment for i i^j and 1 < i, j < k*. 



- n' 




A'2 ■ 




A'm-N-1 


A'm-n 


A'm-2 


A' M-1 




u 


n 






Am-N-2 


Am-n-1 ' 


Am-3 


Am-2 




Vi 




n 




Am-N-3 


Am-n-2 ' 


Am-i 


Am-3 




Vn-1 


Un-1 


Un-2 ■ 


n 




A2 


■ Aff_fj_2 


Am-n-1 




0 




Un-1 ■ 




n 


Ai 


■ ■'^M-Af-S 


Am-N-2 


(18) 


0 


0 


f/„ ■ 


U2 


Ui 


n 




Am-N-3 




0 


0 


0 




Un-1 


Un-2 ■ 


n 


Ai 




. 0 


0 


0 


0 


Un 


Un-1 ■ 




n 





Qa 





Ai 


Am-N-2 


Am-n-1 


Am-n + Un 


■ Am-2 + U2 


Am-1 + U^ 


Am-1 + Ui 


n+A^ 


Am-N-3 


Am-N-2 


Am-n-1 


■ Am-3 + ^3 


Am-2 + U2 


Am-2 + U2 


Am-1 + Ui ■ 


Am-n-^ 


Am-N-3 


Am-N-2 


• + U3 


Am-3 + U3 


Am-N+2 + f«-2 










■ Am-n + Un 


Am-n+1 + Un 


Am-n+1 + Un-1 










Am-n-1 


Am-n + Un 


Am-n + Un 




O+Am 


Ai 


A2 


Am-N-2 


Am-n-1 


Am-n-1 


Am-n + ■ 


■ Am-1 + f 1 


n +Am 


Ai 


Am-N-3 


Am-N-2 


Am-N-2 


Am-n-1 


• Am-2 + U2 


Am-1 + Ui 


H + Am 


Am-n-^ 


Am-N-3 


Az 


A3 


■ Am-n + Un 


Am-n+1 + Un-1 


Am-n+2 + Un-2 


D+Am 


Ai 


I A, 


A2 


Am-n-1 


Am-n + Ufj 


Am-n+1 + Un-1 


■ Am-1 + Ui 


H+Am 



The basic generator of the bulk queue which is concerned with only the arrival and service is a matrix of order 
[M'Z'[likik\] given above in (19) where Qa=Ao + + A2 (20) 
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Its probability vector w' gives, w'Qa =0 and w'. e = 1 (21) 
It is well known that a square matrix in which each row (after the first) has the elements of the previous row 
shifted cyclically one place right, is called a circulant matrix. It is very interesting to note that the matrix 2/i is 
a block circulant matrix where each block matrix is rotated one block to the right relative to the preceding block 
partition. In (19), the first block-row of type [Y]l^ik^k\] x[MY]t=ikik\] is, W = {n + A^.A^, A2, 
...,Af^_i^_2, /1m_jv-i, ^m-n + + ^2' + ^i) which gives as the sum of the blocks {H + 

Am) + Ai+ A2+...+Am-n-2 + ^M-N-i + ^M-N + 1]n+---+^m-2 + ^2+ + ^^1= wWch is the matiix 

given by 

Q"i ^1,2 ^1,3 ■'■ ^l,k* 



Q"= 



■^22,1 Q"2 ■^2,3 ■'■ ^2,k* 
^3,1 ^3,2 Q"3 "■■ ^3,k* 



(22) 



-^k*,l ^k*,2 ^k*,3 ■'■ Q"k*- 

where using (5) and (6), Q'\ = ((7; + T\ a^^iSiIk'.') + ( Ik,<S>iSi + S\ ^J) + diag ((Qi);,;) for 1 < i < k*. The 
stationary probability vector of the basic generator given in (19) is required to get the stability condition. 
Consider the vector w = ( n\(pi ® 0^, n'2(P2 ® 02:- • ^'fc*<?'fc* ® 4'k*) where n' = (n\,n' 2, ... ,n\t') is the 
stationary probability vector of the environment, q)i and (pi are the stationary probability vectors of the arrival 
and service PH processes (7; + 7 j «;) and {Si + S i Pi) respectively. It may be noted n' i{(pi 0 4'i)[{{Ti + 
T'iCCi)®li,') + Uk,mi+S' iPi))] =0. This gives <(<p, ® = n' i{Q^)i_i {cpi ® (t>i) 1 = n'i(Q,)i^i 

{cpi ® (pi) for 1 < i < k*. Now the first column of the matrix multiplication of wfl" is 7r'i((?i)i,i<Pi,i0i,i + 

T^'2{.Qi)2,i<Pii<t>ii[i<P2®<t>2)e]+ + TT'fc. ((?i)fc.,i<Pii0ii[(<Pfc. (8)0fc.)]e = 0 since {(pi®(pi)e = 1 and 

n'Qi=0. In a similar manner it can be seen that the first column block of wQ" is tt:' i{Qi)i i(Pi <Sl (pi + 
^'2 ('?i)2,i'?'i ® (Pi[i'P2 ® 02)e] + + ((?i)fc.,i<Pi O 0i[(<Pfc. O 0fc.)]e = 0 and i-th column block is 

T^'l<iQl)l.i<Pi ® <PiK<Pl ® <Pl)e] +n'2 {Ql)2,l9l ® <Pl{{(P2 ® 02)6] + +Tl'i {Qi)i,i(Pi ® <pi + ...+ 

''^'k*{Qi)k*,i9i® 'PiK'Vk*® 4>k*)]e= 0- This shows that w{n + A^) + wA-^+wA2+...+wAm-n-2 + 
wA^f_fj_i + wA[^_fj + wl]fj+...+wA^f_2 + WU2 + vi/yl^_i + wf/i= w Q"=0. So (w, w,...,w) .W= 0 = (w, w, 
....w) W where W is the transpose W. This shows (w,w...w) is the left eigen vector of and the 

corresponding probability vector is \m''m''m' ' ^'^^'"^ ^ given by 

w = ( 7r'i((pi 0 0i), n'2{(p2 0 (P2), , T^'k*i(Pk* ® 0fc.) ) (23) 

Let q)i = {(pi J ) and (pi = ((pi j ) be the stationary probability components of the arrival and service processes. 
Neuts [5], gives the stability condition as, w' Aq e < w' A2 e where w is given by (23). Taking the sum 
cross diagonally in the Aq and A2 matrices, it can be seen using (9) that 

w' Ao e=^ wQ:^^,nAje=^ {YJil=in=inn' i{(pi ® (g) e) ) = 

^ {i:':i=in=inn'i{cpiTl^®<Pie))=^ (S!i'=iSf;in<(,p,r;,j)=^(Sf;i7r',i:«=ini:'ii ^l^j) 
=^{i:ti n't 1.^=, n Y.jU <Pi.j (T'i)j (>„)) =^(n*, n'l t^U <Pi.j ( T' t)j ^qx)<w'A2 e 
4 w(ILinUje=^ {rn=in=inn'iicpi<^cPi)ie ® S,', J ) =i (S=ii:?=in<(,?^,e ® 0,5;,J ) 
4 (ILiIf=in<(0i5;,j)=^(Zf^i<S=inS;ji0y (SlJj) 

=^(I]?=i7r'iI]!y=ini:^''J^0y (S' i)jQqn))=^Q:'l*^n'iY'^^i(l>i,j(S' i)jEqxp). This gives the stability condition 

as Z^:,n\Y.%,cpij{T\)jE(jx) <Y.^=,n\Zj:,cPij{^^ (24) 
From the result (24) result for the case of PH/PH/1 bulk queue without environment Ramshankar et al., [8] can 
be deduced. This result (24) is the stability condition for the random environment PH/PH/1 bulk queue with 
random sizes of arrivals and of services where maximum arrival size is greater than the maximum service size. 
When (24) is satisfied, the stationary distribution of the queue length exists Neuts [5]. Let n(0, i, j) : for 1 < i < 
k* 1 <j <fci; n(0, k, i, j, j') ; for 1 <k<M-l; 1 < i < k*; 1 < j < fc;; 1 <j <fc'i and n(n, k, i, j, j'): for 0<k< 
M- 1; l<i<k*; l<j</Cj; I <}<k'i and n > 1 be the stationary probability vectors of Markov chain X(t) states. 

Let7ro=('t(0,l,l),7r(0,l,2),...7i(0,l,fci),7r(0,2,l),7i(0,2,2), 7i(0,2,fc2), ,7i(0,k*,l),7i(0,k*,2), 7i(0,k*,fefcJ, 

71(0,1, l,l,l),7i(0,l, 1,1,2)... 71(0,1, l,/Ci,/c'O,7i(0,l,2,l,l),<0,l,2,l,2)...7i(0,l,2,/c2,fc'2),<0,l,3,l,l),7i(0,l^ 

7l(0,l,3,fe2,fe'2)...Jl(O,l,k*,l,l),7l(O,l,k*,l,2) 7l(0,l,k*,fefc.,feV.),7l(0,2,l,l,l), 7l(0,2,k*,/Cfc.,/cV.) 

7i(0,M-l,l,l,l),n(0,M-l,l,l,2) 7i(0,M-l,k*,/<:fc.,/<:'fc.))bea vectoroftype MY.-liki + (M - l)I]f=ifei k'l]. 

Let 7r„=(7i(n,0,l,l,l),7i(n,0,l,l,2) 7i(n,0,l,fei, /<:'i),7i(n,0,2,l,l),7i(n,0,2,l,2) 7i(n,0,2,/c2, /c'2),7i(n,0,3,l,l), 

7i(n,0,3,l,2)...7t(n,O,3,/<:3,fe'3)...7i(n,O,k*,l,l),7i(n,O,k*,l,2) 7i(n,0,k*,/<:fc.,/<:'fcJ,7i(n,l,l,l,l),7i(n,l, 1,1,2) 

7i(n,l,l,/ci,/c'i),7i(n,l,2,l,l),7i(n,l,2,l,2) 7i(n,l,2,fe2,fe'2),7i(n,l,3,l,l),7i(n,l,3,l,2) 3r(n,l,3,/<:3, fe'3) 

7l(n,l,k^l,l),7I(n,l,k^l,2)...7l(n,l,k^fcfc.,fe'fcJ,7l(n,2,l,l,l),...7I(n,2,k^fefc.,feV.),JI(n,3,l,l,l)... 

International organization of Scientific Research 6 I P a g e 



Ph/Ph/1 Bulk Arrival and Bulk Service Queue With Randomly Varying Environment 



n(n,3,k*/(:fc.,feV.)----Jt(n,M-l, 1,1, l),n(n,M-l, 1,1,2) n{n,M-'^X*,kk* k*)) be a vector of type 

lx[MXf=i fej fe'j]- The stationary probability vector n = (719,71^,713,...) satisfies the equations 



nQA=Q and ne=l. (25) 

From (25), it can be seen nQBi+ 71^82=0. (26) 

n^B^+n^A^+n2A2 = 0 (27) 

A2 = 0, for n > 2. (28) 

Introducing the rate matrix R as the minimal non-negative solution of the non-linear matrix equation 

Ao+RAi+R'^A2=0, (29) 

it can be proved (Neuts [5]) that 7r„ satisfies the following. 7r„ = for n > 2. (30) 

Using (26), ttq satisfies Uq = 71^52 (31) 

So using (27) and (31) and (30) the vector can be calculated up to multiplicative constant since satisfies 

the equation 71^ [S2(-Si)"iSo + ^1 + ^^2] =0- (32) 

Using (31) and (30) it can be seen that 71^ [B2 (-Bi)"^eH-(l-R)-ie] =1. (33) 



Replacing the first column of the matrix multiplier of Tii in equation (32), by the column vector multiplier of 
TTi in (33), a matrix which is invertible may be obtained. The first row of the inverse of that same matrix is Ui 
and this gives along with (31) and (30) all the stationary probabilities of the system. The matrix R is iterated 
starting with R{0) = 0; and finding R(n + l)=-AoAj^-R^(n)A2Aj'^, for n > 0. The iteration may be 
terminated to get a solution of R at a norm level where ||i?(n -I- 1) — /?(n) || < e. 

2.3. Performance Measures of the System 

(i) The probability of the queue length S = r > 0, P(S=r) can be seen as follows. For 1 < r < M-1, P(S =r) = 
Xf=i Ii^j=i ^^2=1 '■'ji'j^y For r > M, let n and k be non negative integers such that r = n M -H k. Then 
P(S=r) =Xf;i "^("' ^' iJiJi) , where r = n M -h k, n > 1 and k > 0. (34) 

(ii) The probability that the queue length is zero is P(S =0) = I]?=i l!^Li "^(0- ^J')- (35) 

(iii) The expected queue level E(S), can be calculated. Using (35) and (34), it may be seen that 
E(S)=i:? rP(S = r)=i:f;i HjU 0^(0, ij) +YI^=^ YJtU tj[=i t^U kn{Q, k, 1,^2) 

+I.n=i n=i I-jU ^'2=1 ^'A.;2)(nM + fc) 

=I:L"i' T.i=i t'jl^i kniO, k, t,A,;2) + I.n=i 7r„.(Mn,. . .,Mn,Mn+l,. . .,Mn+l,Mn+2,. . .,Mn+2,. . ., 
Mn+M-l,...,Mn+M-l) = SLl' ^ Sf^^i S^^^i 7r(0, fe, t,A J2) + MK=i n7r„ e+7ri( / - i^^^^ 

Here ^0, ... 0,1, ...,1,2, ...,2, ...,M - 1, ...,M - 1)' is of type [(Sf^i /Ci/c'i)M]xl column vector 
in which consecutively (2]f=ifeife'i) times 0,1,2,3.., M-1 appear. Let it be called ^' when 0 appears 
(Ei'ili ki)timss and others in that order appear (I]f=ifeife'i) times. 
Then E(S)= 7ro^'-i- 7ri( / - -I- M7ri(/ - i? )-2e (36) 

(iv) Variance of S can be derived. Let r] be column vector r|=[0, . . ,0, 1^, ... 1^ 2^,.., T? , ... (M — 1)^, ... , (M — 
1)2' oi type \{^i=lk*kik'i)Vi\Y.V\x\ which consecutively (i=lk*kik'i)i\vc&i, squares of 0,1,2,3.., M-1 appear. Let 
it be called r|' when 0 appears (Si=\ ^i) times and others in the same manner as in r] appear ( Xf=i fei^'i) times. 
Then it can be seen that the second moment, 

E(52)=i:S' r2 P(5 = r)=Yl[U t^U 0^(0' W') +^f=i ^i=i P,=i ^j2=i ^' Uh.h)k^ 

+1^=1 SLV n(n, k, i.AJiXnM + kf =7roii' + E^^^ n(n - l)7r„ e + Y.n=i nn^ e] + 

I»=i7r„r7 +2MK=in7r„f. 

So, '^{S^^=nQv;+M\n^{l - Ry^2R e ^ n^{l - Ry^e\ + n^{l - Ry^-q ^2Mn^{l - Ry^i, (37) 
VAR(S)=E(52) - [^(^il^may be written from (36) and(37). 

III. MODEL (B) MAXIMUM ARRIVAL SIZE M < MAXIMUM SERVICE SIZE N 

The dual case of Model (A), namely the case, M < N is treated here. (When M =N both models are 
applicable and one can use any one of them.) The assumption (vii) of Model (A) is changed and all its other 
assumptions are retained. 

3.1.Assumption 

(vii) The maximum arrival size M= maX]^<;<|(fj maxj^ <;ij. is less than the maximum service size 
N=maxi<i<fc. maxi <,■ <fc.. ]N . 
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3.2.Analysis 

Since this model is dual, the analysis is similar to that of Model (A). The differences are noted below. The state 
space of the chain is as follows presented in a similar way. 

The state of the system of the continuous time Markov chain X (t) under consideration is presented as follows. 
X(t) = {(0, i, j) : for 1 < i < k* 1 < j < fc,)} U {(0, k, i, j, j') ; for 1 < k <N-1; 1 < i < k*; 1 < j < fc- 1 < j < fc'J 
U {(n, k, i, j, j'): for 0 < k < N-1; 1 < i < k*; \ < ] < kc, I < ] < k\ and n > 1}. (38) 
The chain is in the state (0, i, j) when the number of customers in the queue is 0, the environment state is i for 1 
< i < k*and the arrival phase is j for 1 < j <fej. The chain is in the state (0, k, i, j, j') when the number of 
customers is k for 1 < k < N-1, the environment state is i for 1 < i < k*, the arrival phase is j for 1 < j < fe; and 
the service phase is j' for 1 < j' < fe'j. The chain is in the state (n, k, i, j, j') when the number of customers in the 
queue is n N + k, for 0 < k < N- 1 and 1 < n < oo, the environment state is i for 1 < i < k*, the arrival phase is j for 
1 < j <fcj and the service phase is j' for 1 <j' < fe'j. When the number of customers waiting in the system is r, 
then r is identified with (n, k) where r on division by N gives n as the quotient and k as the remainder. The 
infinitesimal generator Qg of the model has the same structure given in (4) but the inner matrices are of different 
orders. 



Qb= 



0 
0 
0 



B 



A2 

0 
0 



0 

A 2 
0 



0 
0 

Ai 
A 2 



0 

Ao 



(39) 



In (39) the states of the matrices are listed lexicographically as 0, 1, 2, 3, .... For partition purpose the zero states 
in the first two sets of (38) are combined. The vector 0 is of type 1 x [I]f=i fe; + (iV — 1) Xf=i fe; ^'J and is 

0=((0,l,l),(0,l,2),...(0,l,fei),(0,2,l),(0,2,2),...(0,2,fe2), (0,kM),(0,k*,2),. . .(0,k*,fefcj,(0,l,l,l,l),(0,l, 1,1,2) 

....(0,1,1,1, fe'i),(0,l,l,2,l),(0,l,l,2,2)....(0,l,l,2,fe'i),(0,l,l,3,l)....(0,l,l,3,fe'i) (0,l,l,fei,l)... 

(0,l,l,fei,fe'i),(0,l,2,l, 1),(0,1,2,1,2)... .(0,1,2,1, fe'2),(0,l,2,2,l),(0,l,2,2,2)....(0,l,2,2,fe'2),(0,l,2,3,l)....(0,l,2,3, 
fe'2)....(0,l,2,fe2,l)...(0,l,2,fe2,fe'2),(0,l,3,l,l)...(0,l,3,fe3,fe'3)...(0,l,k*,l,l),....(0,l,k*,fefc.,feV),(0,2,l,l,l), 
(0,2,1,1,2).. .(0,2,k*,fe,„fe',.),(0,3,l,l,l)...(0,3,k*,fe,.fefcJ,(0,4,l, 1,1). ..(0,4,k*,fefc.fe'fc.)...(0,N-l,l, 1,1)... 
(0,N-l,k*,fej.j, fe'/jj)) and the vector n is of type \x[NYJt=ikik\] and is given in a similar manner as follows 
n=(n,0,l,l,l),(n,0,l,l,2)....(n,0,l,l,fe'i),(n,0,l,2,l),...(n,0,l,2,fe'i),...(n,0,k*,l,l),...(n,0,k*,fefc.,feVJ,(n,l,l,l,l) 

....(n,l,k*,fefc.,feVJ,(n,2,l,l,l)....(n,2,k*,fefc.,feV.) (n,N-l,l,l,l),(n,N-l,l,l,2) (n,N-l,k*,fefc., fe'^J). 

The matricesS\ and A\ have negative diagonal elements, they are of orders [I]f=i ki + {N — 1) Xi=i fe; ^'J and 
[iV2f=i fe; fe'J respectively and their off diagonal elements are non -negative. 
The matrices A'^andA^ have nonnegative elements and are of order [ iV Xf=i fe; fe'; ] • The matrices 
B'q and B'2 have non-negative elements and ai-e of types [I]f=i k^ + (N — 1) Xf=i k^ fe'J x [ iVXf=i fe; k\ ] 
and [ NYJtU^ik'i ] x fe; + (iV - 1) Sf^^ fe; fe'J and they are given below. Using Model (A) for 

definitions of Aj andA'j , and Uj,Vj , and U and letting O and Q' as in Model (A), the partitioning matrices are 
defined as follows. The matrix B'qIs same as that of A'q with first zero block row is of order[ Xf=i fc;] x[ 
iVl]f=i fei fe'i]. The matrix B'2 is same as that of A '2 except the first column block is of type [ N Xi=i fej fe'i ] x 
[ I.'[=i kt] and is (f/^, 0, ... .0) ' where 



0 



0 



^'0 



0 




^2 ® ^2,N 






0 


0 




0 


Ik 






0 




0 






0 


- 0 


0 


0 


0 


0 


... 0- 


0 


0 


0 


0 


0 


... 0 


Am 


0 


0 


0 


0 


... 0 


Am-1 


Am 


0 


0 


0 


... 0 


A2 


A3 


- A„ 


0 


0 


... 0 


y A, 


A2 


■■■ Aff_i 


Am 


0 


... 0. 



4. 



^k*,N 



(40) 



(41) 
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A', 



A\ = 



B\ = 





1 Vn-2 


■■■ ^3 


U2 












0 U„ 


Un-1 


■■• f/4 


U3 


U2 










u u 




TJ 


I J 


I J 

"3 










0 0 


0 


■■■ Ue 


Us 


u. 










0 0 


0 


II., 


II., , 


II., „ 










0 0 


0 


0 


"N 












. 0 0 


0 


... 0 


0 












n 


Ai 


A2 




Am 


0 


0 


0 


0 ■ 


^1 


n 






A:,t 1 
-'*M— 1 




0 


0 


0 






n 






Am-1 


Am 


0 


0 


J J 

'Jn-m-1 


'JN-M-2 


TJ 

UN-M-3 




n 


A 


A 

^2 


A 

•''m-1 


A 

^'m 




[T 

"N-M-1 


TJ 

"N-M-2 




TJ 


n 


A 


/I 

•''m-2 


•''m-1 


Un-M+1 










Ui 


n 


■■■ •'^M-3 


•'^M-2 


II., o 
'-'N-2 


II., „ 


II., . 
UN-4 




Un-M-2 


Un-M-3 


"N-M- 




/I, 

•"l 




IT 

"N-2 


IT 

Un-3 




IT 

"N-M-1 


IT 

"N-M-2 


IT 

"N-M- 


1 "1 


n 


- n' 


A'l 


A'2 




A'm 


0 


0 


0 


0 ■ 


u 


fi 


A, 




'*M— 1 


/lu 


0 


0 


0 


Vi 


Ui 


n 








Am 


0 


0 


^N-M-2 


Un-M-2 


IT 

Un-M-3 




n 


^1 


A2 


A 

•''m-1 


•''m 


^N-M-1 


^N-M-l 


Un-M-2 






n 


Ai 


■■■ •'^M-Z 


Am-1 




^N-M 


^N-M-l 








n 


■■■ •'^M-3 


Am-2 


Vn-3 




TT 




Un-M-2 


Un-M-3 


Un-m- 


n 

-2 ■■• " 


A 

A\ 




Un-2 


Un-3 




Un-M-1 


Un-M-2 


Un-m- 


-1 ■■■ fi 


12 . 




- Am 


-1 + 


Am + t'Af-M 


Un-M-1 






n + Upj 


- Am 


-2 + Upi-M+r 




Am + ^Af-M 




-M-2 


Un-M-3 




n + Uf^ 


^1 + f/w-i 


^2 + yW-2 




-M-1 


Un-M-2 






/2 + f/ 


N 


Ai + t/«-i 


■■■ ^^ 




Un-M-1 




U2 






/2 + f/„ 


■■• ^^ 



(42) 



(43) 



(44) 



^2 
^3 



^2 



f/, 
^M 
A2 + f^«- 



^3 + 

/I2 + [/„- 



"W-M-1 
Am + f/(Vf_i 



^M + fw- 


M 


M-1 + ^Af- 


M+1 


M-2 + Ufi- 


M+2 


n + Uf^ 









Un-M-1 

Am + Un- 

Am-1 + f Af- 



Ai + U„_ 



The basic generator which is concerned with only the arrival and service is 2g = A'^^ + A \ + A'2- This is also 
block circulant. Using similar arguments given for Model (A) it can be seen that its probability vector is w' = 

(ppp where w is given by w = ( n\((pi (g) ^i), n'2i(P2 ® (^2), , T^'k*i'Pk* ® 0fc.) ) and the 

stability condition remains the same. Following the arguments given for Model (A), one can find the stationary 
probability vector for Model (B) also in matrix geometric form. All performance measures including expectation 
of customers waiting for service and its variance for Model (B) have the form as in Model (A) except M is 
replaced by N. 

IV. NUMERICAL ILLUSTRATIONS 

Numerical cases are presented here to illustrate the application of the study. Three examples are studied namely 
(i) M=N=3 (ii) M=3, N=2 and (iii) M=2, N=3. The environment has two states governed by the Markov chain 

with infinitesimal generator | ^ . In state 1, the arrival time distribution is PH 1 (exponential) with 

parameter 5 and the service time distribution is PH 2with representation ^ ^ ''^j and starting probability 

vector (.6, .4). In environment 2 the arrival time has PH 2 distribution with representation \ ^ "'" j and 

t 2 — 4J 

starting probability vector (.4, .6) and the service time distribution is PH 1 (exponential) with parameter 8. For 
the case (i) M=N=3, in the environment 1, the probabilities of bulk arrivals are taken as \pi=.6, \p2=-2 and 
}p3=.2 and the probabilities of bulk services in the two phases are considered as iq=.5, \q2=-3, }q3=.2, 2Q'i=-4, 
2f72=-6 and 2Q'3=0- In the environment 2, the bulk arrival probabilities are assumed as \pi=A, iP2=-6, iP3=0, 
2Pi=.5, 2P2=-5 and 2P3=0 and bulk service probabilities are taken as 2<7i=.5, 2^2=-^ and 2Q'3=0. For the case (ii) 
M=3, N =2 only two probabilities of case (i) alone are changed namely, \q2=-^ and 1(73=0, keeping the same 
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values of case (i) for all others. For the case (iii) M=2, N=3, the probabilities of case (i) are assumed changing 
only two values namely, iP2=-4 and }p3 = 0. The rate matrix R is taken for calculation after 30 iterations for the 
three models. The results obtained are presented in the table 1 below. The probabilities and expected values 
show significant variation for higher and lower values of M and N. The figm-es (1) and (2) show the variations 
of probabilities. 



Table 1. Results Obtained for the Three Cases. 
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Figure 1. Probabilities of Queue Lengths 
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Figure2.The Probabilities of Customer Blocks 
V. CONCLUSION 

Two PH/PH/1 bulk arrival and bulk service queues with random environment have been studied by identifying the 
maximum of the arrival and service sizes and grouping the customers as members of blocks of such maximum sizes. Matrix 
geometric results have been obtained by partitioning the infinitesimal generator by grouping of customers , environment state 
and PH phases together. The basic system generators of the queues are block circulant matrices which are explicitly 
presenting the stability condition in standard forms. Numerical results for bulk queue models are presented and discussed. 
Effects of variation of rates on expected queue length and on probabilities of queue lengths are exhibited. The decrease in 
arrival rates (so also increase in service rates) makes the convergence of R matrix faster which can be seen in the decrease of 
norm values. The standard deviations also decrease. The PH/PH/1 queue with bulk arrival and bulk service with random 
environment has number of applications. The PH distributions include Exponential, Erlang, Hyper Exponential, and Coxian 
distributions as special cases and the PH distribution is also a best approximation for a general distribution. Further the 
PH/PH/1 queue is a most general form almost equivalent to G/G/1 queue. The bulk arrival models because they have non 
zero elements or blocks above the super diagonals in infinitesimal generators, they require for studies the decomposition 
methods with which queue length probabilities of the system are written in a recursive manner. Their applications are much 
limited compared to matrix geometric results. From the results obtained here, provided the maximum anival and service 
sizes are not infinity, the most general model of the PFI/PFI/1 bulk arrivals and bulk services queue with random environment 
admits matrix geometric solution. Further studies with block circulant basic generator system may produce interesting and 
useful results in inventory theory and finite storage models like dam theory. It is also noticed here that once the maximum 
arrival or service size increases, the order of the rate matrix increases proportionally. However the matrix geometric structure 
is retained and rates of convergence is not much affected. Models with multiple servers with PH distributions may be 
focused for further study which may produce more general results. 
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